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The interface-induced magnetization damping of thin ferromagnetic films in contact with normal- 
metal layers is calculated from first principles for clean and disordered Fe/Au and Co/Cu interfaces. 
Interference effects arising from coherent scattering turn out to be very small, consistent with a very 
small magnetic coherence length. Because the mixing conductances which govern the spin transfer 
are to a good approximation real- valued, the spin pumping can be described by an increased Gilbert 
damping factor but an unmodified gyromagnetic ratio. The results also confirm that the spin-current 
induced magnetization torque is an interface effect. 
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I. INTRODUCTION 

The local magnetization dynamics in a bulk ferromag- 
net are usually well described by a phenomenological 
model formulated in terms of three parameters: Hoff, 
an effective magnetic field; 7, a gyromagnetic ratio; and 
a, a Gilbert-damping constant. The field HeS is a sum 
of contributions from externally applied fields, crystal 
anisotropy, shape-dependent dipolar interactions, and 
exchange interactions which govern ferromagnetic spin- 
wave spectral characteristics. 7 is the ratio of the to- 
tal magnetic moment and the angular momentum of the 
electrons in the ferromagnet; in 3c? transition- metal fer- 
romagnets, such as Fe and Co, it is close to the free- 
electron value 7 « The Gilbert-damping con- 
stant a parametrizes the viscous damping of an excited 
magnetization to the (locally) lowest-energy configura- 
tion. Its value differs considerably for various materials 
and also depends on the temperature and on the impu- 
rity/defect composition of a given sample. The motion of 
the magnetization-direction unit vector m is determined 
by the phenomenological Landau-Lifshitz-Gilbert (LLG) 
equation*! 



dm dm 

— —7m X Hoff -f am X 

dt ' dt 
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The magnetization dynamics of small monodomain fer- 
romagnets are well described by the LLG equation ^ 
down to the micron scale. New effects may play a role 
on the submicron scale, however. The magnetization dy- 
namics is no longer a highly coherent process because in- 
terface and surface roughness are relatively more impor- 
tant in small samples. Many-magnon processes can then 
acquire a sizeable spectral weighlj^ and are observable as, 
e.g. an increased line width of the ferromagnetic reso- 
nance (FMR) ? Another source of additional FMR broad- 
ening is non-local, depending on the environment into 
which the ferromagnet is embedded: a time-dependent 



ferromagnetic order parameter pumps spin currents that 
carry angular momentum (and energy) into adjacent con- 
ducting materials!^ This angular-momentum loss, in 
turn, is equivalent to an additional damping torque on 
the magnetization.® 

The spin-pumping concept for the magnetization dy- 
namics of nanostructures has far-reaching consequences. 
It gives rise to an enhanced Gilbert damping of mag- 
netic films in contact with conducting media^ may be 
employed as an FMR-operated spin battery,-^ explains a 
dynamic exchange coupling in magnetic bilayersp^ as 
well as a dynamic stiffness against current-induced mag- 
netization reversalpifi The analysis of experimental FMR 
probes of the magnetization dynamics in single films^ 
and magnetic bilayerg® relied on phenomenological mod- 
els of the electronic structure. Here we show how these 
assumptions can be relaxed by using instead scattering 
matrices calculated from first principles which take into 
account the detailed atomic and electronic structure of 
the materials under study. 

An early phenomenological treatment of the non- 
locality of the magnetization dynamics in hybrid normal- 
metal/ferromagnet {N/F) structures was given by Sils- 
bee et ali^ Recently, Simanek et al}^ pointed out that 
time-dependent linear-response theory could be used to 
calculate the spin flows generated by a ferromagnet with 
a time-varying magnetization in contact with a non- 
magnetic conductor, as an alternative to the scattering- 
theory approach of Tserkovnyak et al~- In spite of the 
different starting point, complete agreement between the 
two methods was demonstratedi^ for the simple case of 
a (5-function magnetic layer embedded in a free-electron 
gas. In addition, it was argued in Refs. 12 13 that the 
electron-electron interactions can considerably enhance 
the spin currents into normal metals with large Stoner- 
cnhancement factors. The linear-response framework has 
also been used to calculate the enhanced Gilbert damping 
of finite-thickness ferromagnetic films.'ii It was argued 
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there that ultra-thin films display oscillatory damping (as 
a function of thickness) due to quantum-size effects. In 
the following, we show that these quantum-interference 
effects are greatly overestimated by the ballistic free- 
electron band model and do not survive when realistic 
transition-metal band structures are used. By calculat- 
ing from first principles the scattering matrix entering the 
spin-pumping theoryji we show that quantum-size oscil- 
lations are much smaller than those reported in Ref. |l4|, 
especially if even small amounts of disorder are intro- 
duced. We also find that the additional term in the fer- 
romagnetic equation of motion is of the Gilbert-damping 
form, with only a very small correction to the gyromag- 
netic ratio (the same conclusion can also be drawn from 
previous work-'-^). Furthermore, the electron-electron in- 
teraction effects discussed by Simaneki^ are taken into 
account in the exchange-correlation potential which we 
calculate self-consistently within the local spin density 
approximation (LSDA) of density-functional theory. Fi- 
nally, our results confirm that the spin-current-induced 
magnetization torquei^ is an interface effect, which was 
earlier taken for grantediii^ and analyzed in detail in 
Ref.lii 

This article is organized as follows. The general theory 
of spin pumping and its consequences for the dynamics 
of the precessing ferromagnet are reviewed in Sec. ^ In 
Sec, mil we describe the first-principles methods used to 
obtain the results presented and discussed in Sec. IIVI A 
comparison with results based on a free-electron model 
is made in an Appendix and conclusions are drawn in 
Sec. El 



II. THEORY 

We first consider a ferromagnetic film of thickness d 
connected to two perfect non-magnetic reservoirs by two 
leads which support well-defined scattering states. The 
electrons incident on the ferromagnet from a lead are 
distributed according to the Fcrmi-Dirac statistics of the 
respective reservoir, whereas the probability that an elec- 
tron leaving the ferromagnet returns there with finite 
spin (or phase) memory is vanishingly small. Such per- 
fect spin sinks can be realized experimentally by attach- 
ing leads to the ferromagnetic film in the form of point 
contacts with dimensions smaller than the electron mean 
free pathJfi Alternatively, a normal conductor with a 
very high spin-flip to momentum scattering-rate ratio (as 
could be provided by heavy impurities with large spin- 
orbit interaction in a light metal or a heavy metal with 
phonon or defect scatterers) can serve as a good spin 
sinkp^ 

Coherent motion of the magnetization, whose direction 
is given by the unit vector m(i), leads to the emission of 
a spin current 



Is = — RcA'^m X h liaA^^ 

47r \ dt dt 



(2) 



per unit area of the contact into each normal-metal 
layer (4 which we will here assume is then fully ab- 
sorbed by the spin sinks (reservoirs)^ The complex spin- 
pumping conductance^! 



(3) 



is the difference between the reflection {g"^^) and trans- 
mission (5*1) mixing conductances (per unit contact 
area) which are defined in terms of the spin-dependent 
reflection and transmission matrices of the ferromagnetic 
film asi^ 



mn 



mn mn J ' 



(4) 
(5) 



Here, S is the F/N contact area, m and n denote scatter- 
ing states at the Fermi energy of the normal-metal leads. 
For spin-pumping into one of the normal-metal layers, g^^ 
is expressed in terms of the amplitude for an incom- 
ing electron in state m of the normal metal to be reflected 
at the interface with the magnetic film into the outgoing 
state 71, while g^^^ is expressed in terms of the amplitude 

for an incoming electron from the other N layer to 
be transmitted across the ferromagnet into the outgoing 
state n. The total angular- momentum loss of the ferro- 
magnet is given by a sum of contributions Q from the 
two leads, characterized by two spin-pumping parame- 
ters and As explained in Ref. 0, adding this 
source of spin angular-momentum current to the right- 
hand side of Eq. (QJ leads to a new LLG equation for the 
monodomain thin film with saturation magnetization Mg 
embedded in the non-magnetic conducting medium, with 
the modified constants aeff and 7cff 



1 _ 1 
7cff 7 
Toff 
7 
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(6) 
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It can be easily shown^ that the real part of A^^ is always 
non-negative so that the correction to the damping is 
always positive. The reader is referred to Sec. IIVI for 
a discussion of the absolute and relative values of g^^ 
and (7*1 . Anticipating these results, we note here that 
in typical situations g*^^ and Img^^ (and thus Imyl^-'') are 
negligible so that the only effect of the spin pumping is to 
make an additional contribution to the Gilbert-damping 
parameter. We shall therefore assume for the rest of the 
current section that <C g^^ with the latter quantity 
being essentially an interface property. 

Eq. ||2Jl was derived for an N/ F/ N structure connected 
to perfectly equilibrated reservoirsi^ By using this ge- 
ometry, the finiteness of the Sharvin conductances is au- 
tomatically included>2^ To apply calculated mixing con- 
ductances to the discussion of spin transport in diffuse 
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systems which are not ideal spin sinks, the "bare" con- 
ductance has to be corrected^i for the correspond- 
ing "spurious" Sharvin resistance as discussed in Ref. 
Additionally, a non-vanishing backflow and reabsorption 
of the spins emitted by the ferromagnet has to be taken 
into account. The latter can be achieved by considering 
the diffusion equation for the spin accumulation in the 
normal lead with Eq. |(2Jl providing the boundary condi- 
tion (see Ref. 0). This leads to an effective (complex) 
conductance A^^ (for either interface) entering equations 
© and |7J) where 



1 



1 

9n 



2e^ 
h 



RsD 



tanh(L/AsD) ' 



(8) 



is the Sharvin conductance of the normal-metal layer, 
given by the number of the transverse channels per spin 
and unit area of the interface;-'- -Rsd = Asd /o" is the unit- 
area resistance of the normal-metal film with conductiv- 
ity cr/2 (per spin) and thickness Asd, the spin-diffusion 
length; L is the actual thickness of the normal-metal 
layer. The last term on the right-hand side of Eq. © 
accounts for impurity, defect, or phonon scattering in the 
normal metal. (Scattering in the ferromagnet on length 
scales longer than the transverse spin-coherence length 
does not modify the result.) When spin- flip scattering in 
the N layer vanishes, Asd oo, A^^ (z.e the back- 
flow spin-current completely cancels the pumping effect) 
and the magnetization dynamics is not modified at all. 

A similar analysis can be applied to magnetic damp- 
ing in more complex multilayer systemsi^- For example, 
in an F/N/F structure the presence of two ferromag- 
netic layers can make damping possible for each individ- 
ual layer even in the absence of spin-flip relaxation in 
the system. In this case, each magnetic layer acts as the 
sink for the spin current pumped by the other layer. If 
the structure is weakly excited from a coUinear equilib- 
rium state, and the individual ferromagnetic resonances 
are well separated, then a different effective conductance 
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enters Eqs. JHl and |7J. Instead of the sum 



A 



for the two magnetic films, the quantity ^lf,„r/r with 



^F/N/F 



2e2 
h 



F/N/F 

L 



(9) 



defined for the globally diffuse system should be used, 
where g^^ is the mixing conductance for the i-th F/N 
interface. Eq. @ can be intuitively interpreted in terms 
of resistances in series: in order to be absorbed, the spin 
current must be pumped through the first F/N interface 
((7|j^ renormalized by 2gf^), propagate across the nor- 
mal layer (L/a term) and enter the second ferromagnet 
through the other interface {g^^^ renormalized by 2g^^ 



N 



The formula for A^p/j^/p can be straightforwardly de- 
rived using the spin-diffusion approach of Ref. |^ It is 
worthwhile pointing out that it remains correct for non- 
diffusive normal metal spacers (cr ^ oo) if the interface 
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mF (int-l) 
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2.742 
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1.410 


mjv(int-l-l) 
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0.036 


mF(int-l-l) 


2.687 


2.691 


1.545 


1.540 


rriF (int-|-2) 


2.336 


2.396 


1.635 


1.596 


mF (int+3) 


2.325 


2.363 


1.621 


1.627 


rriF (iiit-|-4) 


2.238 


2.282 


1.627 


1.624 


rriF (bulk) 


2.210 


2.210 


1.622 


1.622 



TABLE I: Layer-resolved magnetic moments in Bohr magne- 
tons for single N/F interfaces (A'^^Au, Cu; _F=Fe, Co). 



disorder is sufficient to suppress any quantum-size effects 
(see Ref. |25(). 

The effect of spin-dependent scattering on the time 
evolution of the magnetic order parameter is therefore 
mostly governed by three parameters: the reflection and 
transmission mixing conductances of the ferromagnetic 
layer, g"^^ and g^^, and the Sharvin conductance of the 
normal metal, gf^. We noted before^ that these quan- 
tities are in principle accessible to ab initio electronic- 
structure calculationsii^iSSiSi In the following we demon- 
strate this by studying two representative N/F material 
combinations: Au/Fe(001) and Cu/Co(lll), the former 
routinely used by the Simon-Fraser group22i22i2£ and the 
latter by the Cornell group . 



III. FIRST-PRINCIPLES METHOD 

Parameter-free calculations of transmission and reflec- 
tion coefficients were performed using the local spin den- 
sity approximation (LSDA) of density-functional theory 
(DFT) in a two-step procedure. In the first step, the self- 
consistent electronic structure (spin densities and poten- 
tials) of the system was determined using the layer TB- 
LMTO (tight-binding linear mufSn-tin orbital) surface 
Green's function (SGF) method in the atomic-sphere ap- 
proximation (ASA).'^^ The exchange-correlation poten- 
tial in the Perdew-Zunger^^ parametrization was used. 
The atomic-sphere (AS) potentials of 4 monolayers on 
either side of the magnetic layer (or interface) were it- 
erated to self-consistency while the potentials of more 
distant layers were held fixed at their bulk values. Be- 
cause both of the systems we consider, Au/Fe(001) and 
Cu/Co(lll), are nearly ideally lattice matched, common 
lattice constants were assumed for both metals of a given 
structure: ocu/Co = 3.549 A and aau/Fo = \/2 x 2.866 
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=4.053 A. In the second step, the AS potentials serve as 
inputs to calculate scattering coefficients using a recently- 
developed scheme based on TB-MTOsi2S'2Zi2£ Disorder is 
modeled by allowing a number of interface layers to con- 
sist of NxFi-x alloy which is modeled using repeated 
lateral supercells. Because a minimal basis set of s, p 
and d orbitals is used, we are able to treat lateral su- 
percells containing as many as 200 atoms in which the 
two types of atoms are distributed at random in the ap- 
propriate concentration. For disordered interfaces, the 
AS potentials were calculated self-consistently using the 
layer CPA approximation in which each layer can have a 
different alloy compositioniS 

Little is known from experiment about the atomic 
structure of metallic interfaces. We model "dirty" in- 
terfaces with one (for N/F/N systems) or two (for single 
N/ F interfaces) atomic layers of a 50 %-50 % alloy. Such 
a model is probably reasonable for Cufcc/Cofcc because 
of the nearly perfect lattice match and structural com- 
patibility. The situation is however more complicated for 
Aufcc/Fcbcc because of the large difference in AS sizes for 
Au and Fe with Wigner-Seitz radii of 2.99 and 2.67 Bohr 
atomic units, respectively. We have assumed here that 
the disorder is only substitutional and that the diffused 
atoms occupy the AS of the same size as that of the host 
element. In the Au/Fe/Au case, where the alloy is only 
1 atomic monolayer (ML) thick, we assume that the Fe 
atoms diffuse into Au. While the validity of this model 
can be questioned, the insensitivity of the final results to 
the details of the disorder (e.g. one versus two monolay- 
ers of alloy) indicate that this is not a critical issue. The 
layer-resolved magnetic moments for single interfaces are 
given in Table ^ They agree well with values reported 
previously in the literaturei^SiSiM 

The two-dimensional Brillouin zone (2D BZ) summa- 
tion required to calculate the mixing conductances using 
Eqs. Q and lO was performed using A;||-mesh densities 
corresponding to 10^ points in the 2D BZ of a 1x1 inter- 
face unit cell. The uncertainties resulting from this BZ 
summation and from impurity ensemble averaging are of 
the order of a few times 10^^ 17~^m~^, which is smaller 
than the size of the symbols used in the figures. 
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FIG. 1: Reflection spin-mixing conductance (per unit area) of 
a Au/Fe/Au(001) trilayer with perfect interfaces as a function 
of the thickness d of the Fe layer. In this and subsequent 
plots, mixing conductances expressed in terms of number of 
conduction channels per unit area are converted to 0~^m~^ 
using the conductance quantum e^/h, i.e. G^i = {e^/h)gii. 
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FIG. 2: Transmission spin-mixing conductance of a 
Au/Fe/Au (001) trilayer with perfect interfaces as a function 
of the thickness d of the Fe layer. 



IV. RESULTS AND DISCUSSION 

Figures □ toll show how Gf^ = (eV/i)5tx and = 
(e^//i)(7*l depend on the thickness d of the magnetic layer 

(measured in atomic layers) for specular (fc|| -preserving) 
Au/Fe/Au(001) and Cu/Co/Cu(lll) systems. Both 
quantities exhibit oscillatory behavior with, however, no- 
ticeably different periods and amplitudes. The values of 
both G|j^ and G*| are determined by two factors: the 
matching of the normal metal and ferromagnetic metal 
states at the interface (described by the scattering coeffi- 
cients of the single interface) and the phases accumulated 
by electrons on their passage through the magnetic layer 



(quantum-size effect). The first factor determines the 
amplitudes of the oscillations and (for G^^) the asymp- 
totic values, while the second is responsible for the ob- 
served periodicity. In order to better understand this, it 
is instructive to interpret the transmission and refiection 
coefficients of the finite-size magnetic layer in terms of 
multiple scattering at the interfaces. We first note that 
both Cu and Au have only one left- and one right-going 
state at the Fermi level for each value of fc|| and spin so 
that the summations in Eqs. Q and lO reduce to in- 
tegrations over the 2D BZ involving the complex- valued 
functions ?'°'(fc||) and t°'{k^). Retaining only lowest-order 
thickness-dependent terms, dropping explicit reference to 
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FIG. 3: Reflection spin-mixing conductance of a Cu/Co/Cu 
(111) trilayer with perfect interfaces as a function of the thick- 
ness d of the Co layer. 
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FIG. 4: Transmission spin-mixing conductance of a 
Cu/Co/Cu (111) trilayer with perfect interfaces as a function 
of the thickness d of the Co layer. 



fc|| and to the primes on t , we then have 



f 



\cr,a 



where 



(10) 
(11) 



p = [ti, . . . ,t'^) is a vector of transmission 
coefficients between a single propagating state in the nor- 
mal metal and a set of states in the ferromagnet, A'^ is a 
diagonal matrix of phase factors e^^i^"^ {j is an index of 
the states in the ferromagnet), r^^^ is a scalar reflec- 
tion coefficient for states incoming from the normal metal 
and Tp^p is a square matrix describing reflection on the 
ferromagnetic side. The set of states in the ferromagnet 
consists of both propagating and evanescent states. The 
contribution of the latter decreases exponentially with 
the thickness of the layer. 

Concentrating first on the thickness dependence of 
we notice that, in view of Eq. (|HJ|I . the summa- 
tion in Eq. (jSJ is carried out over terms containing 



phase factors e*^'^»^~'^j^'''^. Because of the large differ- 
ences between majority and minority Fermi surfaces of 
the ferromagnet, this typically leads to rapidly oscillat- 
ing terms which mostly cancel out on summing over k\ | . 
It can be arguedi^ in the spirit of the theory of inter- 
layer exchange coupling^S that the only long-range con- 
tributions originate from the vicinity of points for which 
~ ^j±) ~ 0' corresponding to the stationary 
phase of the summand in Eq. These contributions 
will then exhibit damped oscillations around zero value 
as seen in Figs. |5| and 0] 

Turning to g^^, we find on substituting Eq. Ijlll) into 
Eq. Q that there are two thickness-independent contri- 
butions. The first comes from summing the Snm term 
in Eq. and is nothing other than the number of 
states in the normal metal (i.e. the Sharvin conduc- 
tance). The second comes from the rj^^jv'^iv^TV term 
and provides an interface-specific correction to the first. 
Superimposed on these two is the contribution from the 
thickness-dependent terms which, to lowest order, con- 
tain phase factors e^'-'^^^'^'^i^^'^ and e~''^''''^+''iJ-^'^. Just 
as in the case of g^^, one can argue that the integral over 
these terms will have oscillatory character. However, the 
oscillations will have different periods and occur around 
the constant value set by the first two contributions. It 
is clear that the value approached asymptotically by g^^ 
is simply the reflection mixing conductance evaluated for 
a single interface. 

The period and damping of oscillations of g^^ and g*^^ 
as a function of the magnetic-layer thickness d clearly de- 
pends (through the A"') on the electronic structure of the 
internal part of the magnetic layer, which for metallic sys- 
tems is practically identical to that of the bulk material. 
The amplitudes, on the other hand, are related to the 
interfacial scattering coefficients introduced in Eqs. H1U|I 
and Hll|) . Analyzing the scattering properties of the sin- 
gle interface enables us in the following to understand 
why the amplitudes of oscillation of g*^^ are substantially 
larger than those of g^^ for the two systems considered. 
We begin by noting that the transmission probability 
for states in the majority-spin channel assumes values 
close to one over large areas of the Brillouin zone for 
both Cu/Co and Au/Fe, as illustrated in Fig. |SJa) for 
the Cu/Co(lll) interface. For Cu/Co, this results from 
the close similarity of the corresponding Cu and Co elec- 
tronic structures. The situation is more complicated for 
Au/Fe because the majority-spin Fermi surface of Fe con- 
sists of several sheets, unlike that of Au. However, one of 
these sheets is made up of states which match well with 
the states in Au. In the minority-spin channel, on the 
other hand, the transmission probability varies between 
and 1; see Fig. \^h). The maximum size of the (abso- 
lute value of the) "spin-mixing" products of Eqs. Q and 
Q are therefore determined mostly by the majority-spin 
scattering coefficients while the modulation, as a function 
of , is governed by the corresponding minority-spin co- 
efficients. 
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FIG. 5: Plotted within the first Brillouin zone for the Cu/Co(lll) interface: transmission probability for (a) majority- and 
(b) minority-spins, (c) Real and (d) imaginary parts of rj^^^rj^.^jy . (e) Real and (f) imaginary parts of tj^^t^'int where 
^int = i-F^N ■ ^N^F discussed in the text. Note the different scales for panels (a), (b) and for (c) - (f). 

The small reflectivity for the majority-spin states has 
a direct consequence for the values of the mixing con- 
ductances. In the case of g^^ , the second term under the 
sum in Eq. Q will typically have a negligible magnitude. 
This follows directly from r]^_,j^ ~ and Eq. and 

is illustrated in Fig. |S1 (c) and (d) for the r]v^^^rjy'_^^ 
term. As we can see, the only non-zero contributions in 
this case come from the outer regions of the Brillouin 
zone, where states from the normal metal are perfectly 
reflected because of the absence of propagating majority- 
spin states in the ferromagnet. Independently varying 
phases (as a function of for "up" and "down" reflec- 
tion coefficients lead, in the course of integration over , 
to additional cancellation of already small contributions. 
The final outcome is that the values of g^^ are determined 
mostly by the first term in the Eq. i.e. the Sharvin 
conductance of the lead. 

Because the interface transmission in the majority-spin 
channel is uniformly large almost everywhere in the Bril- 
louin zone, the transmission through the magnetic layer 
also remains large for arbitrary thicknesses and its mag- 
nitude (but not its phase) is only weakly modulated by 



the multiple scattering within the layer. The magnitude 
of the tH-^* product is then modulated mostly by the 
variation of the transmission in the minority-spin chan- 
nel, as a function of fc|| . To demonstrate the effect of the 
interface scattering on y^^, values of the product tl^-^^tl*^ 
are shown in Figs.EJe) and (f) for a Cu/Co (111) inter- 
face. Here, tj^^ is defined as the scalar product of the 
interface transmission vectors: i^^^ = tp^j^ ■ t^f^p. As 
one can see, the values assumed by the real and imagi- 
nary part of this product vary strongly throughout the 
Brillouin zone. Unlike the case of g"^^^, however, the val- 
ues span the entire range from -1 to -1-1. An imbalance 
of positive and negative contributions is therefore more 
likely to produce a sizeable integrated value. The com- 
plex values of t^t'^* are further modified by thickness- 
and fc| I -dependent phase factors discussed above, which 
leads to the oscillatory damping seen in Figs. |2] and 0] 
We compare the magnitude and damping of these oscil- 
lations with those derived from a free-electron model in 
an Appendix. 

Figures and [T] show the same quantities (G|| and 
G'i,) calculated in the presence of disorder modeled by 1 
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FIG. 6: Spin-mixing conductances of a Au/Fe/Au (001) tri- 
layer with disordered interfaces as a function of the thickness 
d of the Fe layer. 
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FIG. 7; Spin-mixing conductances of a Cu/Co/Cu (111) tri- 
layer with disordered interfaces as a function of the thickness 
d of the Co layer. 
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Au/Fe 


clean 


0.40 


0.08 


0.466 


0.005 


0.46 


0.83 


0.46 


(001) 


alloy 


0.39 


0.18 


0.462 


0.003 








Cu/Co 


clean 


0.42 


0.38 


0.546 


0.015 


0.58 


0.46 


1.08 


(111) 


alloy 


0.42 


0.33 


0.564 


-0.042 









TABLE II: Interface conductances in units of 10^^ Q ^m 



ture. The results of such calculations are listed in Tablellll 
for clean and disordered interfaces. The disorder here was 
modeled by 2 ML of 50 % alloy. In spite of this difference, 
the values are practically identical to the asymptotic ones 
seen in Figs. nJElEl and[7| In particular, ImG^^ assumes 
values two orders of magnitude smaller than ReG||, with 
the latter being close to the Sharvin conductance of the 
normal metal. This approximate equality results once 
again from a combination of amplitude (small |r^|) and 
uncorrelated spin-up and spin-down phase effects. 

The values given in Table ^differ somewhat from ones 
reported previously in Ref. 15. There are two reasons for 
this. Firstly, the calculations in Ref. were performed 
using energy-independent muffin-tin orbitals linearized 
about the centers of gravity of the occupied conduction 
states. The current implementation-LSS uses energy- 
dependent, (non-hnearized) MTO's, calculated exactly 
at the Fermi energy which improves the accuracy of the 
method. Secondly, on performing the 2D-BZ integration 
in Eq. it was assumed in Ref. llSl that the contribution 
to the sum of /c|| points for which there are no propagat- 
ing states in the ferromagnct should be neglected. How- 
ever, the lack of propagating states in the ferromagnet 
does not necessarily prohibit the transfer of spin angular 
momentum which can be mediated by evanescent states, 
for example in the case of a magnetic insulator. The con- 
tribution from such points should be included in the 
2D-BZ integration. 



monolayer of 50 % alloy added on each side of the mag- 
netic layer. For both systems we have used 10 x 10 lateral 
supercells. The thickness d in this case is that of the clean 
ferromagnetic layer. For both material systems, the ef- 
fect of disorder is to strongly reduce the amplitudes of the 
oscillations. The reflection mixing conductance becomes 
practically constant at the level of its asymptotic (i.e. in- 
terfacial) value. For G*|, the oscillations are not entirely 
damped out but their amplitude is substantially reduced. 
In fact, the values of G*^ become negligible compared to 
ReG|| for all but the thinnest magnetic layers. In ad- 
dition, we expect that diffusive scattering in the bulk of 
the magnetic layer, which for simplicity has not been in- 
cluded here, will have a similar effect. 

In view of the above results, we conclude that in a 
typical situation A^^ « g^^, where can be calculated 
simply for an N / F interface instead of a complete struc- 



Comparison with experiment 

In Ref. l28l Urban et al. reported room-temperature 
(RT) observations of increased Gilbert damping for a sys- 
tem consisting of two Fe layers separated by a Au spacer 
layer. The magnetization of the thinner of the two ferro- 
magnetic layers precesses in the external magnetic field. 
The other ferromagnetic layer, with the direction of its 
magnetization fixed, acts as a spin sink. No modification 
of the damping coefficient was measured for configura- 
tions without a second Fe layer. The latter finding is 
consistent with the prediction given by Eq. (jSJl in the 
AsD oo limit (well fulfilled for Au) as discussed in 

secnn 

In the presence of a second Fe layer, Eq. should 
be used. Neglecting ImA^^.j^.p leads to jcS — 7 and the 
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FIG. 8: Enhancement of the Gilbert damping coefRcient for 
an Fe/Au/Fe trilayer as a function of 1 jd where d is the thick- 
ness of the excited Fe layer. The filled circles (•) are the RT 
values measured in Ref. and the open one (o) is a low 
temperature value from Ref. Isfll . The theoretical predictions 
based on Eq. (EU for K (with a — > oo) are shown as solid 
and the RT-corrected (with phonon scattering) ones as dashed 
lines. The results of K calculations for a Au/Fe/vacuum 
system are given by crosses (x) and stars (*) for specular and 
disordered interfaces, respectively. The value of the Gilbert 
damping a for a single Fe film is marked with an arrow. 



damping enhancement 



inM.d 



(12) 



where a ~ 0.0046 is the damping measured for a single 
layer. Using^° 7 = 2.1^b/^ and the values of the inter- 
face and Sharvin conductances from Tabled (we assume 
that the values are the same for both Au/Fe interfaces), 
Eq. I|12|l is compared with the experimental dat£^ in 
Fig.|Slfor various assumptions about ct in l^. In the low 
temperature limit and neglecting the residual resistivity 
of the Au layer, a 00, Eq. H12|) yields the solid line 
which is seen to overestimate the damping enhancement 
compared to the measured results. Using finite values of 
a will lead to lower values of A^^ and indeed, it was found 
experimentally^^ that lowering the temperature (increas- 
ing the conductivity) increases the damping by as much 
as about 20 % (open circle in Fig.lS)). If we use the room 
temperature (RT) conductivity due to phonon scattering 
in crystaUine bulk Au,''^ Cph = 0.45 x 10^ fJ"im~\ the 
dashed line is obtained which, as expected, is closer to 
the RT measurements. The agreement with experiment 
can be further improved by taking into account the pos- 
sibility of non-negligible residual resistance*'^ of the Au 
spacer. Assuming, for example, CTics — 0.45x10^ ft^^m^^ 
and CTph = would obviously yield the dashed line in the 
figure while taking Crcs — fph — 0.45 x 10^ il^^m~^ and 
1/a = l/cres + l/cph wiU yield a line very close to the 
measured points. 



The theoretical results represented by the straight lines 
in Fig. |S|are based upon the asymptotic, single-interface 
value of G|| from Table ITU assuming G|j^ to be zero. To 
study possible size-dependent corrections in thin films, 
the experimental system needs to be represented by a 
more realistic model than the symmetric N/ F/ N struc- 
tures discussed in the previous section. The Au/Fe/GaAs 
structure used in Ref. 28 differs from these in two im- 
portant respects. First, the transmission mixing con- 
ductance (5*1) is identically zero because of the insulat- 
ing substrate. Secondly, because the reflection is perfect 
for both spin channels, the thickness-dependent terms in 
Eq. Hll() have larger amplitudes, leading to more pro- 
nounced oscillations of g^^^ than those seen in Figs. Q] 
and 131 To estimate the variation which can result from 
size-dependent corrections, we have performed a series 
of calculations for a Au/Fe/vacuum structure, using vac- 
uum instead of GaAs for simplicity. The mixing conduc- 
tance for the other Au/Fe interface is kept at its asymp- 
totic value f Table Hll) . The results for perfect (specular) 
structures;^ marked in Fig. |S1 with black crosses (x), ex- 
hibit oscillations of non-negligible amplitude about the 
asymptotic values given by the solid line (arbitrarily tak- 
ing the low-temperature regime, i.e. cr —> 00 for refer- 
ence) . The introduction of interface disorder (two ML of 
50 %-50 % alloy) yields values for the damping [stars (*) 
in Fig. |S] essentially averaged back to the limit given by 
the single-interface calculations of Table Hll 

We have thus demonstrated that direct first-principles 
calculations can produce values of the damping coeffi- 
cient in the same range as those measured experimen- 
tally. What is more, by taking into account various other 
sources of scattering in the Au spacer and / or quantum- 
size effects, the calculations can be brought into very 
close agreement with experiment. A more definitive 
quantitative comparison with experiment would require 
a detailed knowledge of the microscopic structure of the 
experimental system which is currently not available. 



Spin current induced torque 



The mixing conductances calculated above, which de- 
scribe how a spin current flows through the system in 
response to an externally applied spin accumulation ji 
(defined as a vector with length equal to half of the spin- 
splitting of the chemical potentials |/x| = (/i| — /i|)/2), 
also describe the spin torque exerted on the moment of 
the magnetic layer (see e.g. Refs. ITsllTTlirS 19.22). Con- 
sider for example the situation where the spin accumula- 
tion has been induced by some means in the left lead only 
and the ferromagnet is magnetized along the z axis. The 
spin current incident on the interface is proportional to 
the number of incoming channels in the lead I?^ = -^gfl^/j, 
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whereas the transmitted spin current is given bjH* 



1 



Rcgi, Imgl, 



Img*^^ Reg 
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H (13) 



(gT+gi)/2 

and the reflected spin current by 
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9lt 



Reg'. 



2^ I 



n 



Sh 







gH' 



where g"^ = J2r 
Biittiker conductances. 




^nm l^nmP ^'''^ the Conventional Landauer- 
Therefore g*^^ determines the 
transverse component of the transmitted spin current 



subject to precession and absorption within the mag- 
netic layer. Similarly, the real and imaginary parts of 
dH" - 9u = J2mn rlnnr^n are related to the components 
of the reflected transverse spin current. The rapid decay 
of (?*| (and g^ — g^^) discussed in previous paragraphs 
as a function of increasing magnetic layer thickness im- 
plies that the absorption of the transverse component of 
the spin current occurs within a few monolayers of the 
N/ F interface. In particular we find that the presence 
of the disorder improves the effectiveness of the absorp- 
tion. The limit g*\-^ — > and g'!^^ g^ corresponds to 
the situation where all of the incoming transverse polar- 
ized spin current is absorbed in the magnetic layer. The 
torque is then proportional to the Sharvin conductance 
of the normal metal. As demonstrated in Figs. 1 1141 El 
and[7| this is the situation for all but the thinnest (few 
monolayers) and cleanest magnetic layers. 



Material Dependence 

The input parameters of spin pumping theory are scat- 
tering matrix elements which are computed using the ef- 
fective potential of Kohn-Sham theory. This potential is 
calculated self-consistently and includes electron-electron 
interaction effects via an exchange-correlation potential 
approximated using the local spin density approximation, 
and the Hartree potential. In particular, the modifica- 
tion of interface parameters as a result of magnetic mo- 
ments being induced in the normal metal by proximity 





Cu 


Ta 


Pd 


Pt 


D(eF)(states/Ry-atom-spin) 


2 


10 


15 


12 


[1 ~ D{eF)hc]-^ 


1.1 


1.9 


4.4 


2.2 




0.56 


0.97 


0.62 


0.68 



TABLE III: Density of states at the Fermi level, Stoner en- 
hanement factor and typical Sharvin conductances for bulk 
fee Cu, Pd and Pt and bcc Ta. Typical values of the Stoner 
parameter, I^c, were taken from references lisi and i46i. 



to a ferromagnet (discussed in the Appendix of Ref ill 
is already included in our results in a self-consistent and 
non-perturbative manner (see Table P). For the Cu and 
Au normal metals we have considered, this effect is small. 
Expressed in terms of a Stoner enhancement, this is re- 
lated to the low Fermi level densities of states of these 
metals, D{ep). Viewing it in this way poses the question 
of the possibility of finding much larger effects for ma- 
terials such as Pd and Pt which have a large density of 
states at the Fermi level (see Table 1111)1 and are known 
'^o be close to a ferromagnetic transition as expressed by 
the susceptibility enhancement x/xo = ['^ — D{£p)Ixc\~^ , 
also included in the Table. To calculate this factor, typi- 
cal values of the Stoner parameter, I^c: were taken from 
references and 14^ 

To examine whether enhancements of the Gilbert 
damping parameter recently reported^2ii£ for thin lay- 
ers of Ta, Pd and Pt compared to Cu are related to their 
large Fermi level densities of states, we need to reexamine 
how the electronic structure enters our description of the 
Gilbert damping. In the spin-pumping formulation, the 
quantities determining the damping enhancement are not 
densities of states but transmission and reffection mix- 
ing conductances determined from the scattering matrix. 
These, we have seen, can be approximated very well by 
ReG|| which is very close to the Sharvin conductance of 
the normal metal. Values of this quantity are given for 
Cu, Ta, Pd and Pt in the last row of the Table. It is seen 
that the Sharvin conductance changes less than D{eF)- 
More significantly, with a maximum for Ta, the trend 
does not correspond to that observed experimentallyi^Si^ 
Cu ^ Ta — > Pd — > Pt. We believe that the explanation 
should be sought elsewhere, possibly in the increasing 
spin-orbit interaction which will lead to the heavier ma- 
terials behaving as more efficient spin sinks. To examine 
this suggestion in detail requires the reformulation of the 
spin-pumping theory to include spin-orbit interaction. 



V. CONCLUSIONS 

In summary, we have calculated the transmission and 
reflection mixing conductances that govern the non-local 
effects in the ferromagnetic magnetization dynamics for 
two commonly used N/F combinations: Au/Fe and 
Cu/Co. In both cases, the transmission mixing conduc- 
tance is much smaller than the reflection mixing con- 
ductance except for the thinnest magnetic fllms, only 
a few atoms thick. Even for such thin fllms, g*^^ is smaller 
than Furthermore, is more sensitive to disorder, 
even a small amount of which reduces it to zero while hav- 
ing only a small effect on g'^^^^ as shown in Figs. 15171 For 
all thicknesses, Reg^^ 3> liag"^^ and Re(7|^ is very close 
to its interfacial value (i.e. the mixing conductance of 
the inflnitely-thick magnetic film). The general formulas 
© and predict that the spin pumping renormalizes 
both the Gilbert damping (a) and the gyromagnetic ra- 
tio (7) of a ferromagnetic film embedded in a conducting 
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FIG. 9: The real part of calculated for a free electron 
model with ef ~ 7 eV (energy measured from the bottom 
of the parabolic conduction band in the normal metal) and 
various choices of the exchange splitting A. The interlayer 
distance is taken to be the same as for Cu/Co(lll) system. 
The results of the first principles calculations (•) from Fig. |1| 
are included for comparison. 

non-magnetic medium. However, in view of the results 
discussed in the previous section, we conclude that, for 
all but the thinnest and cleanest magnetic layers, the 
only effect of the spin pumping is to enhance the Gilbert 
damping. The correction is directly proportional to the 
real part of the reflection mixing conductance and is es- 
sentially an interface property. We also find that oscilla- 
tory effects are averaged out for realistic band structures, 
especially in the presence of disorder. Re^ii (which de- 
termines the damping enhancement of a single ferromag- 
netic film embedded in a perfect spin-sink medium) is 
usually very close to gf^ for intermetallic interfacesi^ii^ 
(being in general bounded by 2f;^^ according to its defini- 
tion, Eq.^. These results also apply to the spin-current 
induced magnetization reversal in intermetallic systems, 
indicating that the "effective field" correction due to the 
imaginary part of the mixing conductance and bulk con- 
tributions to the torque are very small. 
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APPENDIX A: COMPARISON WITH A FREE 
ELECTRON MODEL 



A combination of interfacial and bulk dephasing mech- 
anisms, discussed in Sec. IIVI ensures that in the asymp- 
totic (thick magnetic layer) limit the spin-pumping 
mixing-conductance, A^^, reduces to the reflection mix- 
ing conductance G^^, with the latter quantity assuming 
values which are predominantly real and equal to those 
determined for a single N/F interface. Thinner layers ex- 
hibit oscillatory behaviour which is most pronounced for 
(Figs.[21and21). The amplitude of oscillation however 
is at most 20% of the asymptotic value of and de- 
creases to less than 5% for layers more than lOML thick. 
This fast decay, found even for clean, fully coherent struc- 
tures, contrasts with results reported in Ref.^^for a free- 
electron model. For thin layers, Mills found the damping 
coefficient oscillated with amplitude in the range of 80% 
of the asymptotic value and, for layers several tens of 
MLs thick, it was still of order 10%. This feature of 
the free-electron model is illustrated vividly in Fig. O by 
comparing Re(G*^) for Cu/Co/Cu(lll) from Fig. 4 with 
the corresponding results calculated for free-electrons. In 
our free-electron calculation, the Fermi energy in the non- 
magnetic material was taken to be 7eV in order to obtain 
the correct value for the Sharvin conductance of Cu and 
the effect of changing the exchange splitting A of the fer- 
romagnet was studied. For A = 2, 4, 6 eV, the amplitude 
of oscillation is much larger and the decay is much slower 
than what we find for the more realistic multi-band elec- 
tronic structures. As might be expected, increasing the 
exchange-splitting from 2 to 6 eV leads to a shorter pe- 
riod and more rapid decay of the oscillations. However, 
in order to mimic the parameter- free result, an exchange 
splitting in the range of 10 eV would be needed (not 
shown in the figure). Such a large value cannot be justi- 
fied either on theoretical or experimental grounds. This 
discrepancy illustrates the difficulty of mapping the com- 
plex electronic structure of transition metals onto single 
band models in a meaningful way. Free-electron mod- 
els do not adequately describe the effectiveness of the 
thickness-dependent "bulk" dephasing in the ferromag- 
nct. What is more, they also cannot reproduce the com- 
plex spin- and fc| | -dependence of the interface scattering 
coefficients (illustrated in Fig. |SJ| which results from the 
mismatch of the normal metal electronic structure and 
the quite different majority- and minority-spin electronic 
structures of a ferromagnetic metal. For single band free- 
electron models, the interface scattering coefficients con- 
tain much less structure and consequently this model fails 
to take into account even qualitatively the dephasing ef- 
fect of the interface. 
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